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Abstract. In this article, wc give some existence and smoothness resuhs for the law 
of the solution to a stochastic heat equation driven by a finite dimensional fractional 
Brownian motion with Hurst parameter H > 1/2. Our results rely on recent tools of 
Young integration for convolutional integrals combined with stochastic analysis meth- 
ods for the study of laws of random variables defined on a Wiener space. 
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1. Introduction 



' The definition and resolution of evolution type PDEs driven by general Holder contin- 

Ph ' uous signals have experienced tremendous progresses during the last past years. When 

r-| . the Holder regularity of the driving noise is larger than 1/2, this has been achieved 

\ thanks to Young integrals [H] or fractional integration [TT] techniques. The more deli- 

cate issue of a Holder exponent smaller than 1/2 has to be handled thanks to rough paths 
techniques, either by smart transformations allowing to use limiting arguments [21 IH E] 
or by an adaptation of the rough paths formalism to evolution equations [6l[71[T0]. Alto- 
\ gether, those contributions yield a reasonable definition of rough parabolic PDEs, driven 

at least by a finite dimensional signal. 

With those first results in hand, a natural concern is to get a better understanding 
^ of the processes obtained as solutions to stochastic PDEs driven by rough signals. This 

Q ■ important program includes convergence of numerical schemes (see [5] for a result in this 

direction), ergodic properties and a thorough study of the law of those processes. The 
current article makes a first step towards the last of these items. 

Indeed, we shall consider here a simple case of rough evolution PDE and see what 
r> ' kind of result might be obtained as far as densities of the solution are concerned. More 

I specifically, we focus on the following heat equation on (0, 1) 

Yt = St^+! St-umY)u)dBi, te[0,T], (1) 

where T > is a finite horizon, St stands for the heat semi-group associated with 
Dirichlet boundary conditions, is a smooth enough initial condition, Fi : L^(0, 1) — t- 
L2(0, 1) and B : [0,T] ^ M'^ is a (i-dimensional fractional Brownian motion with Hurst 
parameter H > 1/2. For this equation, we obtain the following results: 
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(1) Existence of a density for the random variable Yt{^) for any t G (0,T] and 
^ G (0, 1), when Fi is a rather general Nemytskii operator (see Theorem 13.81 for 
a precise statement). 

(2) With the additional assumption that the Fj's are regularizing vector fields, we 
obtain that the density of lt(0 is smooth (This will be the content of Theo- 
rem |i]8]). 

To the best of our knowledge, these are the first density results for solutions to nonlinear 
PDEs driven by fractional Brownian motion. Let us point out that we could have 
obtained the same kind of results for a more general class of equations (operator under 
divergence form, general domain D C M", drift term, Gaussian process as driving noise). 
We prefer however to stick to the simple case of the fBm-driven stochastic heat equation 
for the sake of readability and conciseness. 

Our main results will obviously be based on a combination of pathwise estimates for 
integrals driven by rough signals and Malliavin calculus tools. In particular, a major 
part of our effort will be dedicated to the differentiation of the solution to equation ([1]) 
with respect to the driving noise B and to a proper estimate of the derivative. Since 
the equations for derivatives are always of linear type they lead to exponential type 
estimates, which are always a delicate issue. This is where we shall consider some 
regularizing vector fields Fi in ([T]), and proceed to a careful estimation procedure (see 
Section 14. ip . It should also be noticed at this point that the basis of our stochastic 
analysis tools is contained in the celebrated book [12] by D. Nualart, plus the classical 
reference [13] as far as equations driven by fBm are concerned. 

Here is how our article is structured: Section [2] is devoted to recall basic facts on 
both pathwise noisy evolution equations and Malliavin calculus for fractional Brownian 
motion. We differentiate the solution to ([T]) and obtain the existence of the density at 
Section |3l Finally, further estimates on the Malliavin derivative and smoothness of the 
density are derived at Section HI 

Notation. Throughout the paper, we will use the generic notation c to refer to the 
constants that only depend on non-significant parameters. The constants which are to 
play a more specific role in our reasoning will be labelled ci, C2,... 

For any A; G N, we will denote by C^'''(M; M) the space of functions on M which are 
/c-times different iable with bounded derivatives. For any 7 G (0, 1), = C'^{[0,T];M.'^) 
will stand for the set of (ti- dimensional) 7-Holder paths on [0,T]. 

2. Setting 

One of the technical advantages of dealing with the simple case of a stochastic heat 
equation on (0, 1) is a simplification in the functional analysis setting, based on rather 
elementary Fourier series considerations (notice in particular that the considerations 
of [7] can be avoided). We shall first detail this setting, and then recall some basic facts 
on equations driven by noisy signals and fractional Brownian motion. Throughout the 
section, we assume that a (finite) horizon T has been fixed for the equation. 

2.1. Fractional Sobolev spaces. As mentioned above, we are working here with the 
heat equation in the Hilbert space B := L'^{0, 1) with Dirichlet boundary conditions. 
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The Laplace operator A on S can be diagonalized in the orthonormal basis 



CniO •= v^sin(7m,^) {n G N*), with eigenvalues A„ := vr 



2 2 

n . 



We shall denote by the (Fourier) decomposition of any function y ^ B on this 

orthonormal basis. 

Sobolev spaces based on B are then easily characterized by means of Fourier coeffi- 
cients. We label their definition for further use: 

Definition 2.1. For any a > 0, we denote by Ba the fractional Sobolev space of order 
a based on B, defined by 

B^ := l^y e L\OA) : (y")^ < ooj . (2) 

This space is equipped with its natural norm := ||A"y||| = Xl^i 

also set Boo = C(0, 1). 

The above-defined fractional Sobolev spaces enjoy the following classical properties 
(see mm): 



Proposition 2.2. Let Ba,Boo be the Sobolev spaces introduced at Definition \2.1\ Then 
the following hold true: 

• Sobolev inclusions; If a > 1/4, then we have the continuous embedding 

Ba C Soo. (3) 

• Algebra.' If a > 1/4, then Ba is a Banach algebra with respect to pointwise multipli- 
cation, or in other words 

y ■ < MbMIb^- (4) 

• Composition.- If < a < 1/2, (p E Ba and / : M — )■ M belongs to C^'^, then f{(p) G Ba 
and 

\\f{v)\\B^<Cf{l+MBj. (5) 

Here, f{ip) is naturally understood as f{(p){^) := f{(p{^)). 

Let now St be the heat semigroup associated with A, and notice that if an element 
y e L^(0, 1) can be decomposed as y = ^„>i2/"e„, then Sty = ^„>i e~'^"*|/"e„. The 
general theory of fractional powers of operators provides us with sharp estimates for the 
semigroup St (see for instance 



Proposition 2.3. The heat semigroup St satisfies the following properties: 

• Contraction.- For all t > 0, a > 0, St is a contraction operator on Ba- 

• Regularization.- For all t E (0,T], a > 0, St sends B on Ba and 

\\SM\Bc < Ca,Tt~"\\(p\\B- (6) 

• Holder regularity. For all t G (0,T], {p G Ba, 

\\St^ - ^\\b < Ca,Tt"'\\(p\\B„- (7) 
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2.2. Young convolutional integrals. The stochastic integrals involved in equation ([T]) 
will all be understood in the Young sense. In order to define them properly, let us first 
introduce some notation concerning Holder type spaces in time. To begin with, for any 
a > and any subinterval / C [0, T], set Ba) for the space of continuous Sa-valued 
functions on /, equipped with the supremum norm. Then Holder spaces of Sa-valued 
functions can be defined as follows: for k G (0, 1), set 

C^il; B^) ■.= [ye C%I; B^) : sup \~ < oo 

L s<tei r ~ '^1 

Observe now that the definition of our stochastic integrals weighted by the heat semi- 
group will require the introduction of a small variant of those Holder spaces (see [TUl E] 
for further details): we define Ba) as 

C^{I- Ba) := \y e C\I- B^) : sup \\y^-f^~^f\^^ < oo 

In order to avoid confusion, the natural norms on the spaces C^{I\Ba)i C^lI^Ba) are 
respectively denoted by Af[-; C'^{I; Ba)], Af[-] C^{I] Ba)], etc. For the sake of conciseness, 
we shall often write ^^(S^) (resp. instead of C'^([0, T]; (resp. C^([0, T]; 

We also need to introduce a family of spaces iS^) in the following way: 

Lemma 2.4. For any n G (0,1) and any subinterval I C [0,T], let €^''^{1; B,^) be the 
space associated with the norm 

Ml-; C'^'V, B.)] := ^i-, C\I- B,)] + C'^il; B,)]. 

Then the following continuous embedding holds true: 

C°'''(/;i3«) cC''(/;i3). (8) 

More generally, for every X > n, 

Af[y; C^{I- B)] < M[y; C^; B^)] + c. III'-"" Af[y; C\l- By)]. (9) 
Proof. Indeed, owing to ([7]), one has, for every s < t E I, 

\\yt - ysh < \\yt - Stsysh + Wi^ts - ld)2/s||s < \\yt - Stsyshx + cx\t- s\^ Wyshx- 

□ 

With those definitions in hand, the following proposition (borrowed from [^) will be 
invoked in the sequel in order to give a meaning to our stochastic integrals weighted by 
the heat semigroup: 

Proposition 2.5. Consider a '-f-Holder real-valued function x defined on [0,T]. Let 

I = [ii,i2] be a subinterval of [0, T] and fix hi E [0, 7] such that 7 + k > 1. Suppose that 
z G C°(/; B\) n C'*(/; B\_a) for some parameters A > 0, < a < min(K, A). Then, for 
every s < t E I, the convolutional Riemann sum 

converges in B\ as the mesh of the partition Tl := {s = t^ < ti < . . . < tn = t} tends to 
0, and we denote the limit by St-uZudXu. 



MALLIAVIN CALCULUS FOR FRACTIONAL HEAT EQUATION 



5 



Moreover, for every ip G Bx, there exists a unique path y G C'^{I', Bx) such that yn^ = ip 
and yt — St^sUs = jl St-u{zu) dxu if s < t E I. For this function, the following estimate 
holds: 

Af[y-&'{I;Bx)] < c\\x\\,{Af[z;C''{I;Bx)] + Ur'' Af[z;C^{I; Bx^a)]} , (10) 
for some constant c that only depends on (7, n, A, a). 

2.3. Malliavin calculus techniques. This subsection is devoted to present the Malli- 
avin calculus setting which we shall work in, having in mind the differentiability prop- 
erties of the solution to ([T]). 

2.3.1. Wiener space associated to fBm. Let us first be more specific about the probabilis- 
tic setting in which we will work. For some fixed H E (1/2, 1), we consider (f2, J-", P) the 
canonical probability space associated with the fractional Brownian motion with Hurst 
parameter H. That is, Vt = Co([0, T]; M*^) is the Banach space of continuous functions 
vanishing at equipped with the supremum norm, J-" is the Borel sigma-algebra and P is 
the unique probability measure on Q such that the canonical process B = {Bt, t E [0, T]} 
is a d-dimensional fractional Brownian motion with Hurst parameter H, with covariance 
function 

B[BlBi] = 1 {t'^ + - It - s\'^) s,t E [0,T]. (11) 

In particular, the paths of B are almost surely 7- Holder continuous for all 7 E (0, H). 

Consider then a fixed parameter H > 1/2, and let us start by briefly describing the 
abstract Wiener space introduced for Malliavin calculus purposes (for a more general 
and complete description, we refer the reader to [131 Section 3]). 

Let (ei, . . . , erf) be the canonical basis of M*^, S be the set of M'^-valued step functions 
on [0,T] and "H the completion of S with respect to the semi-inner product 

(l[o,t] ei, l[o,s] ej)u ■■= Rh{s, t) l{i=j), s, t E [0, T]. 

Then, one constructs an isometry K*^ : Ti L^([0, T]; M"^) such that Kff{l[o^t] Cj) = l[o,t] 
Knit, ■) Ci, where the kernel K = Kh is given by 

K{t, s) = chs^-" I {u- s)"-iu"-^ du 

J s 

and verifies that E[i?*i?^] = J^^* K{t,r)K(s,r) dr, for some constant ch- Moreover, let 
us observe that K"^ can be represented in the following form: 

[K*^V], = r^rdrK{r,t) dr = dHt""-''^ [n'^-'"^)] , (12) 

Jt ^ -I * 

where stands for the fractional integral of order a. The fractional Cameron-Martin 

space can be introduced in the following way: let KLh '■ -^^([0, T]; M'^) — I-Lh '■= 

1Ch{L'^(\Q,T]]W^)) be the operator defined by 

[lCHh]{t):= [ K{t,s)h{s)ds, h E L\[0,T];R'^). 
Jo 

Then, T-Lh is the Reproducing Kernel Hilbert space associated with the fractional Brow- 
nian motion B. Observe that, in the case of the classical Brownian motion, one has 
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that K{t, s) = l[o,t](s), K* is the identity operator in L'^{[0, T]; M*^) and Tin is the usual 
Cameron-Martin space. 

In order to deduce that {Q, "H, P) defines an abstract Wiener space, we remark that 
H is continuously and densely embedded in Q. In fact, one proves that the operator 
TZh ■ H — )■ Hh given by 

nHiJ:= [ K{-,s)[K*tlj]{s)ds 
Jo 

defines a dense and continuous embedding from "H into f2; this is due to the fact that 
TZnip is if-Holder continuous (for details, see [131 P- 400]). 

Let us also recall that there exists a c?- dimensional Wiener process W defined on 
{fl, "H, P) such that B can be expressed as 

Bt = [ K{t,r)dWr, te[0,T]. (13) 
Jo 

This formula will be referred to as Volterra's representation of fBm. 

2.3.2. Malliavin calculus for B. Let us introduce now the Malliavin derivative operator 
on the Wiener space (fi, P). Namely, we first let S be the family of smooth functionals 
F of the form 

F = f{B{h,),...,B{K)), 

where hi, . . . , hn & H, n > 1, and / is a smooth function having polynomial growth to- 
gether with all its partial derivatives. Then, the Malliavin derivative of such a functional 
F is the "H-valued random variable defined by 

i=l * 

For all p > 1, it is known that the operator V is closable from L^^Q) into Lp{Q; "H) 
(see e.g. [121 Chapter 1]). We will still denote by V the closure of this operator, whose 
domain is usually denoted by ©^'^ and is defined as the completion of S with respect to 
the norm 

\\F\kp ■■= mF\n + mmfn))' ■ 

Sobolev spaces D^'^ for any k E N and p > 1 can be defined in the same way, and 
we denote by Of^^ the set of random variables F for which there exists a sequence 
{Qn,Fn)n>i C T X 3'''^ such that Qn t ^ a-s- and F = Fn a.s. on Qn- We also set 

Remark 2.6. For F G 3^''^, one can write VF = Yl'j=i'^'^ F ej, where V^F denotes the 
Malliavin derivative with respect to the j^^ component of B. 

Since we deal with pathwise equations, we shall also be able to differentiate them in a 
pathwise manner. The relation between almost sure and Malliavin derivatives has been 
established by Kusuoka, and we quote it according to [121 Proposition 4.1.3]. 

Proposition 2.7. A random variable F is said to be 'H-differentiable if for almost all 
u E Q and for any h E H, the map v H- F(ijj ■\- uTZHh) is differentiable. Those random 
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variables belong to the space ^l^^, for any p > 1. Moreover, the following relation holds 
true: 

{VF, h)n = DF{B){nHh), hen, (14) 
where we recall that T) stands for the Malliavin derivative and D for the pathwise differ- 
entiation operator. 

Stochastic analysis techniques are widely used in order study laws of random variables 
defined on a Wiener space. Let us recall the main criterions we shall use in this direction: 

Proposition 2.8. Let F be a real-valued random variable defined on J-", P). Then 

(i) If F e D]'^^ for p > 1 and \\DF\\-}i > almost surely, then the law of F admits a 
density p with respect to Lebesgue measure. 

(ii) If F G anc? £[111)^11^^] is finite for allp > 1, then the density p of F is infinitely 
differentiable. 

3. Existence of the density in the case of Nemytskii-type vector fields 
In this section, we first consider a general equation of the form 

yt = St<f+ [ St.umy)u) dxl , L\0, 1) ,te [0,T], (15) 
Jo 

driven by a rf- dimensional noise x = {x^, . . . ,x'^) considered as a function with 7 G 
(1/2,1). We shall be able to handle the general case of a perturbation involving so-called 
Nemytskii operators: an operator N : B ^ B is said of Nemytskii type if there exists a 
continuous and bounded function / : M — M such that 

iv(^)(o = Njicpm ■■= fi^iO) , ^ G i3 , e G (0, 1). 

In the sequel, we identify the function / with its associated operator Nf. 
Thus, Equation (fT5|) can here be written as 

yt = Stv+ [ St-uif^iyu)) dxi , ifeB,te [o,T], (16) 

Jo 

where the /j : M — M are smooth enough functions. It is readily checked that if each 
fi belongs to C^'*'(M; M) and y G C°'''(-B«) for some k G (max(l - 7, 1/4), 1/2), then the 



integral in the right-hand-side of f lT6|) can be interpreted with Proposition 12.51 Indeed, 
owing to ([5]), we know that f{y) G C^{Bk), while, due to the embedding (|8]), 

Af[fiyy,C^iB)] < \\f'\UAf[y;C^{B)] < c\\f'\UAf[y;C''^{B^)] < 00. 

For the remainder of the section, we shall rely on the following regularity assumptions. 

Hypothesis 1. We consider k, G (max(l — 7, 1/4), 1/2) and an initial condition (p G 
B^. The family of functions {/i, . . . , fd} is such that fi is an element of ^^'''(M; M) for 
i = 1, . . . ,d. 

In this context, the following existence and uniqueness result has been proven in [TJ 
Theorem 3.10]: 

Proposition 3.1. Under HypothesisUl Equation l[T5\) interpreted with Proposition \2. 51 
admits a unique solution y G C^'^^^Bk), where we recall that the space C^''^{Bk) has been 
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As a preliminary step towards Malliavin differentiability of the solution to ([I]) , we shall 
study the dependence on x of the deterministic equation (1T6|) . 

3.1. Differentiability with respect to driving noise. For equation (fT6|) . consider 
the apphcation 

^■.C ^C^'^iB^), x^y, (17) 

for a given initial condition ip. We shall elaborate on the strategy designed in [13] in 
order to differentiate $. Let us start with a lemma on linear equations: 

Lemma 3.2. Suppose that {x,y) G x €^''^[3^) and fix t^ G [0,T]. Then for every 
w G C^''^{[to,T]; Bk) , the equation 

vt=wt+ [ St^uUliVu) ■ Vu) dxi , te [to, T], (18) 

Jto 

admits a unique solution v G C°''^([to, 7"]; -B^), and one has 

Ar[t;;C°'''([to,T];i3,)] < a,,,T ■ Ar[i/;; C°'''([to, T]; i3,)], (19) 

where Cx,y,T '■= C{\\x\\'y,J\f[y;C'^''^{Bi^)],T) for some function C : (]R+)^ — t- M+ growing 
with its arguments. 

Proof. The existence and uniqueness of the solution stem from the same fixed-point 
argument as in the proof of Proposition 13.11 (see [3 Theorem 3.10]), and we only focus 
on the proof of ffT^ . 

Let / = [(ii,(!-2\ be a subinterval of [to^^"]- One has, according to Proposition 12.5^ 
M[v-C^{I;B,)] 

< U[w- r (S.)] + c l/r-'^ ||x||, {U[fl{y) ■ v; C\l; B.)] + Aflf^iy) ■ v; C^{I; B)]} . (20) 
Now, by using successively ([4]) and ([5]), we get 

Af[fl{y)-v;C'iI;B^)] < cM[f:iyy,C%B,)W[v;C'{I;B^)] 

< c{l + ^^[y■,c"{B,)]}^^[v■,c'>iI■,B,)], 

while, owing to dS]) and ([3]), 

Af[f:{y)-v;CV,B)] 

< ^[fM, CV, B)]Ar[v; C\I- B^)] + U[f[{y)-, C\l- B^)]U[v- C^{I; B)] 

< c { 1 + Af[y- C^'^iB,)] } Af[v- C''V, B.)] . 
Going back to fl2U]) . these estimates lead to 

^^[v; cv, B^)] < A^[^; c%B^)] + c,,y ^^[v; c''-{i; b^)], 

and hence 

Af[v; C^'^'il- B^)] < \\ve, + cAf[w; C^(S,)] + c,,, Ur'' Af[v- C^'^'il- B,)]. 
Control f|T9|) is now easily deduced with a standard patching argument. 

□ 

The following lemma on flow-type linear equations will also be technically important 
for our computations below. 
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Lemma 3.3. Fix {x,y) & x and for every u G [0,T], consider the system of 

equations 

^,^ = S,.^{f,{y^))+ f S,.M-iyv.)-K,u)dxi , teKT], zG{l,...,m}. (21) 

J u 

Then, for every i G {1, . . . , m} and t G [0, T], the mapping u ^ ^ is continuous from 
[0,t] to Bk- In particular, for every C, G (0, 1), u i— "^l^i^) is a continuous function on 

Proof. Let us fix z G {1, ... , m}, t G [0, T]. For any < u < v < t, set 

r;j^):=n.-n. , se[v,T]. 
It is easy to check that F* ^ is solution of the equation on [v,T] 



Jv 

Therefore, according to the estimate ( IT^ . 
Now, observe that 

/•V 

J u 

and since y G C^'^'iB^), it becomes clear that ^ 0. 

□ 

We now show how to differentiate a function which is closely related to equation fll6p . 
Lemma 3.4. The application F : x C'^''^{Bk) — )■ C^''^{Bk) defined by 

F{x, y)t := yt - Stip- / St^uifiiVu)) dx^, 



is differentiahle in the Frechet sense and denoting by DiF (resp. D2F) the derivative 
of F with respect to x (resp. y), we obtain 

D,F{x,y){h)t = - f St-u{f^{yu))dhl (22) 

D2F{x,y)iv)t = Vt~ [ St^uUl{yu)-Vu)dxl (23) 

Besides, for any x G C , the mapping D2F{x, $(a;)) is a homeomorphism of C^'^{Bk). 
Proof. One has, for every h G G C'^''^{B^), 
F{x + h,y + v)t - F{x, y)t = 

vt- [ St-uU'liyu) ■ vu) dxl - / St.M'^iyu)) dh^ - [rI{v) + R^h, v)] , (24) 



with ^ ^ ^ 

Rl{v):= [ St^szldxl , 4-.= [ dr [ dr'rf:{y, + rr'vs)-vl 
Jo Jo Jo 
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Rt{v,h):= / St^.zldhl , zl := dr f-{y, + rv^) ■ 
Jo Jo 

and we now have to show that 

Af[R'{v) + R'{h, v)- C''^{B,)] = o + Ar[v; C''^{BS ^'^ 

Observe first that M\R\v);&^''{B^)\ < cM\R]{v);C''{B^)\. Thanks to ^ and (E]), we 
get 

Ikslls^ < c // drdr' + rr't; J < c{l + + IK'sIIbJ 

J J\o,\Y 

Besides, owing to (E]) and ([H]) and setting Mts(r, r') = f'liyt + rr'v^ — f-'ins + rr'vg) we 
end up with 



\\zt-Zs\\B < II ll^tsl'^^'^Olle IK'sIIb^ + cllt^t - t'sllelllt'tllB^ + 

< c{\\yt - VsWb + IK'i - Vs\\b] \\vs\\%^ + c||t;t - WsIIb + IK'slUJ • 

+ U[v-C''^\B,)]U[v-C\B,)] 



The estimate (|T0|) for the Young convolutional integral now provides us with the expected 
control M[R^{v)]C^''^{BtS)\ = 0{Af[v;C°''^{Bf,)f). In the same way, one can show that 
Af[R^{h,v);d^'''{B^)] = 0{\\h\\.^ ■ Af[v;C^^%B^)]), and the differentiability of F is thus 
proved. 

Of course, the two expressions f l22|) and fl23|) for the partial derivatives are now easy to 
derive from flM|) . As for the bijectivity of 1^2-^(2^, ^{x)), it is a consequence of Lemma 13 .21 

□ 

We are now ready to differentiate the application $ defined by fll7p : 

Proposition 3.5. The map $ : — j- C^'^^B^) is differentiable in the Frechet sense. 
Moreover, for every x E C"' and h G C°°, the following representation holds: if t E 
[0,T],eG (0,1), 

D^x){hm= fwUOdhl, (25) 







where G C{[0,t]] B^) is defined through the equation 

= St^u{fm^)u)) + r S,.Uf-iHxU ■ K,u) dxl (26) 

J u 

Proof. Thanks to Lemma 13.4^ the differentiability of $ is a consequence of the implicit 
function theorem, which gives in addition 

L)$(a;) = -D2F{x, $(x))"^ o DiF{x, $(x)), x G C^. 
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In particular, for every x,h E C', z := D^{x){h) is the (unique) solution of the equation 

zt= f St^MHx)u))dK+ f St^MH^)u)-Zu)dxl te[0,T]. (27) 
Jo Jo 

If X G C*" and h G C°°, an application of the Fubini theorem shows (as in the proof of [T3| 

Proposition 4]) that the path Zt := /q ^tuO^^u (which is well-defined thanks to Lemma 

13. 3p is also solution of (!27|) . and this provides us with the identification (1251) . 

□ 

As the reader might expect, one can obtain derivatives of any order for the solution 
when the coefficients of equation (fT6|) are smooth: 

Proposition 3.6. Suppose that fi G C°°'''(M;M) for every i G {!,..., m}. Then the 
function $ : — €^''^{8^) defined by ( 1T7|1 is infinitely differentiable in the Frechet 
sense. Moreover, for every n G N* and every x,hi, . . . ,hn G , the path Zt : = 
D"$(x)(/ii, . . . , hn) satisfies a linear equation of the form 

zt = wt+ [ St^^{f^mx)^)-Zu)dxl , tG[0,T], (28) 
Jo 

where w G C^''^{Bk) only depends on x,hi, . . . , hn- 

Proof. The details of this proof are omitted for the sake of conciseness, since they simply 
mimic the formulae contained in the proof of [T3| Proposition 5]. As an example, let us 
just observe that for x,h,k G , the path Zt := D^^{x){h,k)t is the unique solution 
of ^ with 

wt := r St-uifmx)^) ■ D^x){h)^) dK + f St^^Um^)u) ■ D^{x){k\) dhl 
Jo Jo 

+ / St^M7iHx)u)-D<l>{x){h)^-D<^{x)ik)^)dxl 
Jo 

□ 

3.2. Existence of the density. We will now apply the results of the previous section 
to an evolution equation driven by a fractional Brownian motion B = {B^, . . . , B'^) with 
Hurst parameter H > 1/2. Namely, we fix k G (max(l/4, 1 — 7), 1/2) and an initial 
condition (p E B^- We also assume that fi G C^'^{R; M) for z = 1, . . . , m. We denote by 
Y = $(-B) the solution of 

Yt = StV+ [ St^MYu)) dBi , te [0,T]. (29) 

Notice that since if > 1/2 the paths of B are almost surely 7- Holder continuous with 
Holder exponent greater than 1/2. Thus, equation fl2I?l) can be solved by a direct appli- 
cation of Proposition 13.11 Moreover, one can invoke Proposition 13.51 in order to obtain 
the Malliavin differentiability of Yt{C,): 

Lemma 3.7. For every t G [0,T],^ G (0, 1), Yt{^) G ©j^^ and one has, for any /i G "H, 

{V{Y,{0), h)^ = D^{B){nHhUO- (30) 
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Proof. According to fll4p . we have that 

ae\e=o 

Furthermore, Proposition 13.51 asserts that $ : — > €'^''^{3^^) is differentiable. Therefore 
with 

\R{enHh)t{0\ < M[R{enHh)-C\B^)] < cN[R{enHh)-C\B^)] = o{e), 
and hence j^^^^^^{B + eTZHh)t{C) = D^{x){TZHh)t{C}, which trivially yields both the 
inclusion Yt{C) G dJ^^ and expression (l30l) . 

□ 

With this differentiation result in hand plus some non degeneracy assumptions, we 
now obtain the existence of a density for the random variable Yt{^): 

Theorem 3.8. Suppose that for all A G M, there exists z G {1, . . . , c?} such that fi{\) 7^ 0. 
Then for allt G (0, 1] and^ G (0, 1), the law ofYt{^) is absolutely continuous with respect 
to Lebesgue measure. 

Proof. We apply Proposition 12.81 part (i), and we will thus prove that > 
almost surely. Assume then that ||P(yt(,^))||^ = 0. In this case, owing to fl30|) . 
we have D^{B)(7lHh)t{^) = for every h E Ti. In particular, due to fl25|) . one has 
/o^tit(0'^^M — f*^^ every h G C°°. As u ^tu(0 is known to be continuous, it 
is easily deduced that "ifl^uiO = every u G [0,t] and every i G {!,..., rf}, and 
so = "^If-l^) = fi{Yt{^)) for every i G {1, . . .d}, which contradicts our non- vanishing 
hypothesis. 

□ 

4. Smoothness of the density in the case of regularizing vector fields 

Up to now, we have been able to differentiate the solution to (|T6l) when the coefficients 
are fairly general Nemytskii operators. However, we have only obtained the inclusion 
Yt{^) G ©iqc- Additional problems arise when one tries to prove Yt{^) G D^'^, due to bad 
behavior of linear equations driven by rough signals in terms of moment estimates. This 
is why we shall change our setting here, and consider an equation of the following type 

yt = Stv+ [ St-u{L{f,iyu)))dxl , tG[0,r] , ^ e B, (31) 

where x G C^([0, T]; M.'^) with 7 > 1/2, : M ^ R (? G {1, . . . , rf}) is seen as a Nemytskii 
operator (see the beginning of Section [3]), and L stands for a regularizing linear operator 
of B. Let us be more specific about the assumptions in this section: 

Hypothesis 2. We assume that for every i G {1, . . . ,d}, fi is infinitely differentiable 
with bounded derivatives. Moreover, the operator L : B ^ B is taken of the form 

L(0)(O:= fdvU{^,vmv), 
Jo 
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for some positive regularizing kernel U such that Cu '■= mingg(o,i) drj 1/(^,7]) > 0. In 
particular, we assume that L is continuous from B to B\ for every A > 0. 

In this context, the following existence and uniqueness result holds true: 

Proposition 4.1. Under Hypothesisl^ for any A > 7 and any initial condition if G B\, 
Equation / f57]) interpreted with Proposition \2.5\ admits a unique solution in C"'{B\). 

Proof. As in the proof of Proposition 13.1^ the result can be obtained with a fixed-point 
argument. Observe indeed that if y G C'^{I] B\) (/ := [(.i,(-2\ C [0, 1]) and z is the path 
defined by zi,^ = y/,^, Zt - StsZs = St^u{L{fi{yu))) dx'^ (s < t e I), then, according to 
Proposition 12.51 z G C'^{I] B\) and one has 

Af[z; C\I- B^)] < c\\x\\, {M[LU{y)); C\I- Bl^)] + \I\^ M[LU\y)); C\I- S^)] } • (32) 

Now, owing to the regularizing effect of L, it follows that A/'[L(/(?/)); C°(/; < 
||^||£(B,B;,)||/||oo and 

Ar[L(/(2/));C^(/;S™)] < ||L|U(0,^,)||r|UAr[y;C^(/;S)] <cAr[y;C°'^(/;S,)], 

which, together with (132|) . allows to settle the fixed-point argument. □ 



For the sake of clarity, we henceforth assume that T = 1. The generalization to any 
(fixed) horizon T > easily follows from slight modifications of our estimates. 

Moreover, for some technical reasons that will arise in the proofs of Propositions 14.21 
and 14. 3[ we will focus on the case A = 2 + 7 in the statement of Proposition 14.11 In 
other words, from now on, we fix the initial condition ip in the space -62+7. 

4.1. Estimates on the solution. Under our new setting, let us find an appropriate 
polynomial control on the solution to fl3T|) in terms of x. 

Proposition 4.2. Suppose that y is the solution of l[31\) in C'^{B2+^) with initial condi- 
tion if. Then there exists a constant C^j^l such that 

Af[y;C''{[0,l];B2+,)] < C,,^,^ (1 + ||a:||,) (max (||a;||;/M|^||y/^J)'"' . (33) 
Proof. For any N eW, let us introduce the two sequences 

£k = SN,k ■= ^ _^ ^ , -^0 := , 4+1 = (^k+i '■= + ^N,k- 
The first step of the proof consists in showing that we can pick such that for every k, 

< 1- (34) 
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For the latter control to hold at time (i.e., for k = 0), we must first assume that 

1 /2 

N > ||v^||g2+7' observe that for any k, one has, owing to f llOp . 

< c\\x\\,{^[L{f{y));C\[i,J,^^];BT+^)]+el^[L{f{y));C^ 

< c||x||^||L|U(i5,i5,^^){l+£^Ar[y;C^([4,4+i];^)]} 

< c||x||, {l + elAf[y; C^([4, 4+i]; i32+,)] + e'^'Af[y; C°([4, 4+i]; ^2+,)] } 

< cillxll, {l+el^^[y;C^i[hJk+l];B2+,)]+el'^'\\y,J\B,^^] , 

where we have used Q) to get the third inequality. Consequently, if we take N such that 
2ciA^~^||a;||^ < 1 (i.e. N > (2ci||a;||^)^/^), we retrieve 

Ar[y;C^([4,4+i];S2+7)] < 2cM\^ + el\\yejB,+^ (36) 

and hence 

II2/4+1IIB2+7 < l + (l + 4^^)lb4llB2+7- 
From this estimate, if we assume that £^fc||2/4 ||b2+7 — then 

2 II II 2 2 f II II 7 1 o 2 ^fc+l ^ ~^ ~^ ^) 

and is thus proved by induction. Going back to fl5^ . we get, for every k, 

Ar[y;C^([4,4+i];i32+^)] <2ci||x||^ + l. (37) 

By a standard patching argument, this estimate yields 

^[y;C^{[OA];B,^,)] < {2ci||x||, + 1} 1^^"^ 

where K stands for the smallest integer such that ^1^=0^1^ — ^■ 

Finally, observe that 2 > J2k=o^k = ^k=N \^ ^^i*^ ihx\s one can check that K < 
(e^ — 1)A^ < IN . To achieve the proof of fl33|) . it now suffices to notice that can be 
picked proportional to max(||x||y''', ||v^||b2+^)- ^ 

We now consider a linear equation, which is equivalent to f llSp in our regularized 
context: 

zt = wt+ [ St^u{L{f^{yu)-Zu))dxi , tG[0,l], (38) 
Jo 

where w G C^{B2+'y) and y stands for the solution of ( 13T|1 with initial condition ip G B2+-y 

The existence and uniqueness of a solution for ( 138|) can be proved along the same lines 
as Proposition 14. that is to say via a fixed-point argument. We shall get a suitable 
exponential control for this solution. 

Proposition 4.3. There exists constants Ci,C2 which only depend on f, L and 7 such 
that 

Af[z;C'{[0,l];B2+,)] < C,Af[w;C''^{B2+,)]exp (c^max [m'JI^, . (39) 
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Moreover, if Wt = St'ip for some function ip G -82+^ , there exists an additional constant 
C3 which only depend on f , L and 7 such that 

< C3||^|U.,,max [m'JI^, exp (c^ max [m'JI^, (40) 

Proof. We go back to the notation eN,kAk of proof of Proposition 14. 2[ and set, for 
every c > 0, 

N{c) :=max(||^||^f^^,(2c||a:||,)i/^). 

We have seen in the proof of Proposition 14.21 that there exists a constant ci such that 
for every > N{ci) and every /c, one has simultaneously 

<1 , U[y,C\r,A^,^,]-B,^,)]<2c,\\x\\, + l. (41) 
Suppose that > iV(ci) and set Mw '■= A/'[ty; C'^{B2+'y)\. One has, similarly to ([35]), 

< AC + c||x||, {M[ny) ■ z- c°([£f , Ci]; s'")] + 4,.A/'[r(i/) ■ ^; Ci]; ^'")]} 

< + c||x||,{Ar[.; C°([£f , S)] + .;,,Ar[.; C^([£f , ^f^J; S)] 

< AC + c2||x||,Ar[.; c°([£f , Ci]; ^2H-,)] {1 + 4,.A/'[y; , Ci]; ^)]} 

+C2||x||,4,,Ar[z;C^([£f,£f^J;S2+,)], 

where we have used ([3]) and ([2]) to derive the last inequality. Therefore, if we choose 
> max (A^(ci), (2c2||x||^)^/'''), one has, for any N > N2 and any k, 

Ar[z;C^([£f,d;S2+,)] 

< 2AC + 2c2||a;||,A/-[2;C°([£f ,£fj;i32+,)] {l + e^,,Ar[y; C^([£f , ^)]} . 

Thanks to (Q and (14T1) . we know that 

A/'b;C^([^f,C+i];^)] < c{Ar[y;C->([CC+,];S2+,)] + 4,.lbcb.+.} 

< c{2ci||a;||^ + 2} . 

As a consequence, there exists C3 such that for any > 

ATlz; C^([£f , i^^,]; 02+,)] < 2AC + Cs\\x\\,M[z; C\[i^ , ff+J; {l + . 

(42) 

Then, for any N > N2, 

Af[z-C'^{[e^,i^^,];B2^,)] 

< 2AC+ +C3||x||,e;^,Ar[^;C°([£f,Ci];^2+,)] {l + 

Pick now an integer > N2 such that 

C3A^3'^||X||,{1 + A^3-^||X||,}<1, 
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and we get, for any k, 

k 

so Af[z]C^{[0,l]]B2+^)] < 2^(^3)||m;o||b2+^ + 2^(^^)+2AC, where K{Ns) stands for the 

smallest integer such that '^^Iq^^ £N3,k > 1- As in the proof of Proposition 14.2^ one 
can check that K^N^) < cN^. In order to get it suffices to observe that there 

exists a constant C4 such that any integer A^3 > C4 max(||<y9||g^^^, 11x11^'*') meets the above 
requirements. 

Suppose now that Wt = Stip. In particular, AC = 0. Then we go back to fH21) to 
obtain, thanks to fl39p . 



Af[z;C^{[i^\i^l,];B2^,)] < C^Ms^^^N^ exp (c^ max [M'H^, \\x\\l/^ 
which entails 

Ar[z;C\[0A];B2+,)] < Cil|^||i5,^^iV3^i^(iV3)^-^exp(c2max(||vp||^f H-H^/^ 

< CsU'h^^^N.exp (C^max 

and f HOjl is thus proved. 



□ 



Remark 4.4. For any to £ [0, 1], the proof of Proposition 14.31 can be easily adapted to 
the equation starting at time to 

zt = Wt,to + / St-u{,L{f-{yu) ■ Zu)) dx\ , w,^, e C'^{[to, 1]; , t e [to, 1], 

J to 

and both estimates fl39l) and fl40l) remain of course true in this situation. 



4.2. Smoothness of the density. Let us now go back to the fractional Brownian 
situation 

Yt = Sty^+ [ St^u{L{MY^))) dBi , te[0,l], i32+7, (43) 

where 7 G (|, iJ) is a fixed parameter. We suppose, for the rest of the section, that the 
initial condition (p is fixed in and that Hypothesis [2] is satisfied. We denote by Y 
the solution of ( 143|) in C'^{B2+'y) given by Proposition 14.11 

As in Subsection 13. 2[ we wish to study the law of Yt{^) for t G [0, 1] and ^ G (0, 1). 
Without loss of generality, we focus more exactly on the law of Yi{^), for ^ G (0, 1). 

The first thing to notice here is that the whole reasoning of Section [3] can be trans- 
posed without any difficulty to Equation (143|) . which is more easy to handle due to the 
regularizing effect of L, as we have seen in the previous subsection. Together with the 
estimates fl33p and this observation leads to the following statement: 

Proposition 4.5. For every ^ G (0,1), Yi{C,) G D°° and the law of Yi{C,) is absolutely 
continuous with respect to the Lebesgue measure. 
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Proof. The absolute continuity of the law of Yi{C,) can be obtained by following the lines 
of Section [31 which gives ^ well. Then, like in Proposition 13. 6[ observe that 

n-th (Frechet) derivatives of the flow associated with (143|) satisfy a linear equation 
of the form 

= W,^ + f S,.^ ■ Z:)) dBl t G [0, 1]. 

The explicit expression for {n > 1) can be derived from the formulae contained in 
[131 Proposition 5], and it is easy to realize that due to fl5^ . one has N\W'^; C°''^(i32+^)] G 
L^iyt) for any n and any p. Then, thanks to (l39i) . we deduce that N'[Z'^]C^{B2+'y)] is 
a square-integrable random variable, which allows us to conclude that ^ (see 

[la Lemma 4.1.2]). □ 



The following proposition, which can be seen as an improvement of Lemma 13.31 (in 
this regularized situation), provides us with the key-estimate to prove the smoothness 
of the density: 

Proposition 4.6. For every s G [0, 1] , consider the system of equations 

%^, = St^sW,{Ys))) + I' St-um;{Yu) ■ K,s)) dBi , tG[s,l] , zG{l,...,m}. 

(44) 

Then, for every i G {1, . . . , m} and every t G [0, 1], G C'''([0, t]; -62+^). In particular, 
for any ^ e (0,1), ^ C^[0,t]). 

Moreover, one has the following estimate 

■C^([0,t];S2+,)] < QUMb,^,, II^II,) -exp (cmax [M'^^, PII^^)) , (45) 
for some polynomial expression Q. 



Proof. As in the proof of Lemma 13. 3[ we introduce the path 

r;j5):=^l,.-^U , se[v,l],0<u<v<t, 
and it is readily checked that ^ solves the equation on [v, 1] 



F^,,(5) = 5,_.(vi/;,-M/;j + J Ss-ULUWu,)-K^uH))dBi. 

Therefore, thanks to the estimate fl39l) . we get 



< 



cllK,. - nJl6.+.exp (cmax Pll^^)) • (46) 

Then, by writing 

= L(/,(n)-/,(K))-[^,_. - Id] (L(/,(K)))- r S._.(L(/;(F^)-v^j„J) dB: 

J u 



J 

Wi 
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we deduce that 

II ^ v,v ^ 

+ \\L\\c^B,B2,jB\\,{Af[nY) . ^l^;C%B"^)]+Af[nY) ■ ^l^;C^{B"^)])} 

<c\v-u\''{i + \\B\\,}{i + ^^[Y; cm} {1 + A/'[*U; ^"(^2+,)] + ^f[Ku, c^m} ■ 

Going back to (146|) . the result now easily follows from the embedding C^''^' {B2+'y) C C^'iB) 
and the three controls fl33|) . (1391) and (itoj). 

□ 

Proposition 14.61 implies in particular that the Young integral jl^tuiO^h^ is well- 
defined for every h G C^, t G [0, 1] and ^ G (0, 1). We are thus in a position to apply 
the Fubini-type argument of [131 Propositions 4 and 7] so as to retrieve the following 
convenient expression for the Malliavin derivative: 

Corollary 4.7. For every ^ G (0, 1), the Malliavin derivative ofYi{^) is given by 

Vl{Y,iO) = ^UO , 3G[0,1], zG{l,...,m}, (47) 
where \E'*^ stands for the solution of [JW on [s, 1]. 

Theorem 4.8. Suppose that there exists Aq > such that for every i G {!,..., m} and 
every 77 G M, /^(r/) > Aq. Then, for every ^ G (0, 1), the density ofYi{^) with respect to 
the Lebesgue measure is infinitely differentiable. 



Proof. We shall apply here the criterion stated at Proposition 12.81 item (ii). Notice that 
we already know that ^1(0 G so it remains to show that for every p > 2, there 
exists eo{p) > such that if e < eo{p), then P (^i(O) ||k < e) < . 

To this end, we resort to the following practical estimate, borrowed from [21 Corollary 
4.5]: for every (5 > H — 1/2, there exist a > such that 

P {\\T^{yim\n <e)<V {\\V{Y,m\oo < e") + P {\\V.{Y,{0)h > e-") ■ (48) 

The first term in the right-hand-side of (l48l) is easy to handle. Indeed, owing to the 
expression (ITT]) for the Malliavin derivative of (0 5 



\V.{Y,m\oo > . inf = . inf |L(/,(F0)(OI 

2=1, ...,m 1=1, ...,m 



inf 

i=l,...,m 



rfr/?7(e,r/)/,(Fi(r/)) 



> cc/Ao > 



(remember that U and C[/ have been defined in Hypothesis [2]), so that P(||r'.(Fi((^))||oo < 
e") = for e small enough. 

Then, in order to cope with P (||P,(Yi((^))||^ > e""), one can simply rely on the 
Markov inequality, since, according to fH5l) . 



\mY,{0)h = ll*i,.(OII/3 < c sup Af[^l-X''{[0,l];B2+,)] 

ie{l,...,m} 

< cQ (11^11^2+^, ||5||^) ■ exp (cmax (||v?||s2+^, 11^11^^)) 
which proves that ||P.(Yi(^)) ||^ G L'^{Q) for every q>l. 



□ 
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